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Ram pressure stripping and the formation of cold fronts 
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ABSTRACT 

Chandra and XMM-Newton observations of many clusters reveal sharp discontinuities in the 
surface brightness, which, unlike shocks, have lower gas temperature on the X-ray brighter 
side of the discontinuity. For that reason these features are called "cold fronts". It is believed 
that some cold fronts are formed when a subcluster merges with another cluster and the 
ram pressure of gas flowing outside the subcluster gives the contact discontinuity the char- 
acteristic curved shape. While some edges may not arise directly from mergers (e.g., A496, 
[bupke & White 2003), this paper focuses on those which arise as contact discontinuities be- 
tween a merging subcluster and the ambient cluster gas. We argue that the flow of gas past the 
merging subcluster induces slow motions inside the cloud. These motions transport gas from 
the central parts of the subcluster towards the interface. Since in a typical cluster or group 
(even an isothermal one) the entropy of the gas in the central regions is significantly lower 
than in the outer regions, the transport of the low entropy gas towards the interface and the 
associated adiabatic expansion makes the gas temperature immediately inside the interface 
lower than in any other place in the system, thus enhancing the temperature jump across the 
interface and making the "tip" of the contact discontinuity cool. We illustrate this picture with 
the XMM-Newton gas temperature map of the A3667 cluster 

Key words: methods: numerical, galaxies: general, galaxies: clusters: individual: A3667, 
X-rays: galaxies: clusters 



1 INTRODUCTION 



Cold fronts were discovered as sharp features in the X-ray sur- 
face brightness dist ribution in Chandra observations of the clusters 
A2142 a nd A3 667 iMarkevitch et al.'2000: Vikhli nin et alj2001iJ. 
see also IMarkevitch et al. 2002.') . Similar features have now been 
found in several other clusters iSun et all l2002t iKempner et alj 
l2002t ). A natural assumption that these sharp features were due to 
shocks propagating through the gas was immediately rejected since 
the measured gas temperature was lower on the X-ray brighter side 
of the feature. It was suggested instead that the observed edges were 
formed as a result of the flow of hotter gas around the colder gas 
cloud, the front itself being the contact discontinuity separating hot 
and cold gas (Markevitch et al. 2000; Vikhlinin et al. 2001a). The 
motion of the cloud with respect to the gas on larger scales can be 
understood, for example, as the result of a merger of a subcluster 
with a more massive and hotter cluster. And indeed, features resem- 
bling cold fronts were found in numerical simulations of cluster 
formation (Nagai & Krav tsov 2003; Bialeket al. 20Q2^. 

In this short paper we argue that ablation of the gaseous cloud 
causes characteristic differential motion of the gas inside the sub- 
cluster, which transports the low entropy gas from the subcluster 



core towards the contact discontinuity, thus enhancing the jump in 
temperature and surface brightness across the discontinuity. 



2 ILLUSTRATIVE MODEL 
2.1 Numerical Setup 

In order to study the global dynamics and temperature structure of 
cold clouds, we ran a number of hydrodynamic simulations. Initial 
conditions were chosen to corre spond roughly to those observed 
in A3667 t Vikhlinin et all2001ah . however, we varied all relevant 
parameters within an order of magnitude to confirm that the con- 
clusions we draw from the simulations are robust. 

The simulations were performed using the FLASH hydro code 
dPrvxell etalJi200l]h . This code is an adaptive mesh refinement 
code using a high precision PPM solver, explicit in time and for- 
mally accurate to 2nd order. To verify our results, we also ran 
a batch of simulations using the ZEUS 3D code iClarke et alj 
1994; Stone & Norman 1992a b). We restricted most of the sim- 
ulations to two dimensions, the bulk of the simulations assuming 
axi-symmetry. We ran several test cases in a cartesian semi-infinite 
slab geometry to confirm that the axi-symmetric assumption did not 
corrupt our results. Furthermore, we re-ran our fiducial simulation 
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in 3D (at lower resolution) to confirm that the results are not an 
artifact of enforced symmetries. 

We used an adiabatic equation of state. Both codes used incor- 
porate shock capturing schemes. However, to keep the simulations 
scale free, we did not include the effects of radiative cooling. For 
our fiducial simulation presented in the following section, it turns 
out that cooling of the coolest gas is at most marginally important 
over the simulation time span: For our fiducial simulation (see be- 
low) the minimum cooling time anywhere during the simulation is 
icooi.min ~ 4 Gyrs, compared to a total simulation time of 4 Gyrs. 
While this is clearly a simplification, the effects of cooling will only 
increase the temperature contrast between the low and high entropy 
regions and thus enhance the effects discussed in this paper. 

The initial setup consists of a stationary isothermal atmo- 
sphere^ of cold, dense gas and a wind of uniform hot, light gas 
moving at uniform speed relative to the isothermal atmosphere. The 
isothermal atmosphere is assumed to be in hydrostatic equilibrium, 
following a beta model density profile of the form 



where we chose /3 = 0.5, = 250 kpc, and pc = 3.6 x 
10~^ cm~'^. For simplicity, we assume that the gravitational field 
binding the cold gas is provided predominantly by dark matter. 
Thus, we keep the gravitational potential fixed throughout the sim- 
ulations. The uniform gas is assumed to have an electron density 
Px ~ 4.6 X 10~^ cm~^, with a relative velocity such that its inter- 
nal Mach number is M — v^jc^ = 1.1, while the Mach number 
with respect to the sound speed of the cold cloud is Mx = 1.6. We 
used a temperatu re normalisation of 7. 7 keV for the undisturbed 
external medium iVikhlinin et alj2001al) and 3.2 keV for the undis- 
turbed isothermal cloud, which sets the length and velocity scale of 
the simulation. 

The computational volume was restricted to the range of r G 
[0 Mpc, 2 Mpc] and z G [-2 Mpc, 2 Mpc] . Increasing the physi- 
cal dimensions of the simulated volume by a factor of two did not 
affect the internal structure of the cloud significantly. The initial 
interface is located at the z = —500 kpc surface. Inflow boundary 
conditions were chosen for the z = —2 Mpc boundary (forcing ve- 
locity, pressure, and density to v^^Px, and respectively), outflow 
boundary conditions for the z = 2 Mpc and the r = 2 Mpc bound- 
aries, and reflective boundary conditions along the r = Mpc axis. 

The effective (i.e., maximum) resolution within this box was 
varied between 1024 x 512 and 8192 x 4096 without significant 
impact on the global dynamics. The bulk of the simulations were 
run at an effective resolution of 2048 x 1024 (corresponding to a 
cell size of 1.95 kpc), with a refinement depth of 8, while forcing 
the resolution in the central 500 kpc to be at least 3.9 kpc. For a 
typical size of ~ 500 kpc (as observed in A3667) this implies that 
the cloud is safely resolved with more than 100 resolution elements 
across a cloud radius. 

In order to evaluate the observed temperature distribution, 
we calculated emission weighted temperature maps, using the ob- 
served background temperature and background intensity around 



the cold front in A3667 (due to the limited box size and assumed 
uniform external density we cannot calculate the background emis- 
sion from the host cluster self-consistently). We assumed pure 
bremsstrahlung emission, which is a good approximation above 
T ^ 2 keV and sufficient for our purposes. 

2.2 Simulation results 

The simulations generally show the following temporal evolution: 
the interface between the moving and the stationary gas drives a 
strong shock into the stationary, dense gas, while a weak shock 
is driven into the hot gas (depending on the initial velocity, this 
weak shock develops into a compression wave or a weak bow shock 
upstream). The interface between the compressed/shocked gases is 
a contact discontinuity, which envelopes the cold gas in roughly a 
bullet shape. 

The front of this contact discontinuity is roughly spherical, 
while Kelvin-Helmholtz instabilities and turbulence, induced by 
the tangential velocity discontinuity between cold and hot gas, lead 
to the development of large eddies I Maz zotta et a l. 2002i). 

The shock passage through the cold gas moves the material 
with the lowest entropy and the highest density away from the cen- 
ter of the potential well. After the shock passage, this leads to a 
reversal of the velocity pattern, due to the gravitational field whose 
center is no longer co-spatial with the lowest entropy, highest den- 
sity material. This field accelerates the material against the direc- 
tion of the background flow, transporting it back towards the con- 
tact discontinuity. 

As material is streaming towards the front of the cloud, a 
developing shear-layer due to unstable Kelvin-Helmholtz modes 
drives material along the contact discontinuity backwards, in the 
direction of the background flow. Together with the backflow 
induced by the displacement of the cloud core away from the 
center of the gravitational well, this results in the formation 
of a strong vortex inside the cloud, which transports material 
inside the cloud against the direction of the background flow 
along the symmetry axis (i.e., to the left in Fig. 0, and, along 
the surface, in the direction of the background flow (i.e., to the 
right in Fig. 0. It is this vortex that leads to the transport of 
low entropy, low temperature material towards the front of the 
cloud. We have visualised this transport in an animation viewable at 
http : / / www . mpa-garching . mpg . de/^heinzs / transport . htm^ 
Morphologically very similar structures hav e been seen in previou s 
simulations of ram pressure stripping by iBalsara et alj il994h : 
lMurravetai]jl993h . 

The effect of the internal circulation on the entropy and tem- 
perature structure of the cloud can be seen in Fig. 1. The lowest 
entropy material is transported to the front, while material with 
higher entropy fills the center of the cloud. The rarefaction wave be- 
hind the cloud leads to rather low temperatures on the trailing side 
as well, due to the adiabatic expansion of the material in the rar- 
efaction. However, the emission weighted temperature map shows 
clearly that the dense, high emission measure material at the front 
of the cloud produces the lowest observable projected temperature. 



^ We note that the assumption of isothermality is idealistic - any realistic 
halo will presumably possess a stratified atmosphere with low entropy gas 
at the center. However, for the puipose of our study, isothermality is a con- 
servative assumption, as it decreases the temperature contrast that can be 
achieved and more clearly reveals the nature of the dynamical effects we 
will be discussing. 



2.3 Transport processes and cloud destruction 

In the outskirts of the cloud, the background flow can unbind gas 
from the cloud potential as long as the ram pressure exceeds the 
binding energy of the gas, i.e, as long as > \p{'r) cf. This 

sets the initial size ro of the stripped cloud, after the passage of 



Ram pressure stripping and cold fronts 3 




z[Mpc] ^[Mpc] z[Mpcl 2 IMpcl z[Mpcl 

Figure 1. Simulation results for our fiducial run at time steps 0, 1, 2, 3, and 4 Gyrs from left to right. Top panel: 2D slice of the entropy S = 
log [(p/p^/'') / (kcV cm^ j j (colour scale), overlayed with vector representation of the velocity field (red arrow corresponds to a velocity of 1400 km s~^). 
Second panel: Slice through the temperature distribution, overlayed with density contours. Third panel: Emission weighted, projected temperature (assuming 
uniform background at relative level and temperature measured by Chandra), overlayed with surface brightness contours, assuming a side-on viewing angle 
of 90°. Each box is 1 Mpc on a side. Bottom panel: plot of pressure (blue) and density (red) along the axis. 



the wind driven shock wave (roughly corresponding to the second 
panel in Fig.0: 



7Px 



(2) 



where Cs is the sound speed of the unperturbed, cold, isothermal 
cloud, 7 = 5/3 is the adiabatic index of the gas, and v is the veloc- 
ity of the external wind. We define po as the cloud density at ro. In 
the case of the above simulation, ro ~ 2rc. For large velocities or 
small density contrasts, all the material is unbound and stripping is 
effective even in the absence of dynamical instabilities. This is the 
case if px > |pc . In the simulation shown above and, most 
likely in A3667, this condition is not met JVikhlinin et all2001b^ . 
Thus additional dynamical effects must come into play to destroy 
the cloud. 

The removal of material along the clou d surface is facilitated 
by Kelvin-Helmholtz eddies jNulsedligil . The largest growing 
mode has a wavelength of order Amax ~ ro for ro > r^- In the limit 
of Px <C po (which is achieved both in observed cold fronts and the 
simulation shown above), the growth time of a mode of wavelength 



A is ta = ^ ^J^■ A shear layer of thickness A ~ A/27r set up by 
a fully non-linear mode of wavelength A that transports material 
away from the cloud surface with velocity of order < Ux. Since the 
largest unstable mode has Amax ~ ro, the destruction time of the 
cloud is approximately given by 



Tdost 



few X 



27rro ro 

Amax ^x 



' few X 10 



ro 



(3) 



where we defined the ablative index ( to be calibrated by the simu- 
lations. 

The shear layer induces a vortex inside the cloud, which be- 
comes stronger over time as the width of the shear layer grows as 
longer and longer wavelengths become fully non-linear. By con- 
struction, the time it takes for this vortex to transport the lowest 
entropy material from the core to the front IS "Tturnovcr ^ "Tdcst; 
since this is the time for the shear layer to remove most of the mass 
of the cloud, and thus to turn over the entire cloud. It is therefore 
natural to assume that during the late stages of cloud evolution the 
coldest material will find itself close to the front of the cloud. 

Accordingly, from the simulation shown above, we see that the 
cloud is destroyed at about the same time that the low entropy gas 
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has been transported to the front. The numerical experiment yields 
(, ^ 0.08, which is the order magnitude expected from eq. J3j. 



3 ILLUSTRATIVE CASE - A3667 

To further illustrate the relevance of the mechanism, considered 
above, to the actual structure of the observed cold fronts, we calcu- 
lated the temperatur e structure of A3667 us ing XMM-Newton data. 
The front in A3667 jVikhlinin et all2001 J) is the clearest front de- 
tected by Chandra. 

We used four XMM-Newton observations of A3667 (observa- 
tions IDs 0105260301, 0105260401, 0105260501, 0105260601). 
The description of the data and detailed analysis are given by 
iBriel & et alJ i2003h . Here we provide only the gas temperature 
map for this cluster based on the MOS data. To calculate the pro- 
jected gas te mperature di stribution, we employed the method de- 
scribed in Ch urazov et alJ l 1999). Namely, we fit the observed spec- 
tra at a given location as a linear combination of two template spec- 
tra corresponding to emission from an optically thin plasma (con- 
volved with the MOS energy response) with temperatures 2 and 8 
keV and determine the gas temperature as a function of the relative 
weights of the template spectra. The resulting "temperature map", 
adaptively convolved to have an effective number of counts per 
smoothing window of 6400, an approximate S/N of 80, is shown 
in Fig.|2| The contours show isophotes of the X-ray surface bright- 
ness in the 0.5-3 keV energy band. The outerm ost contour ap- 
proxim ately traces the discontinuity, discussed in IVikhlinin et alJ 
From this figure it is clear that the lowest entropy gas is 
located along the interface and its temperature is lower than that in 
any other place. The extension of the cool region outside the dis- 
continuity (marked by the contour) is due to the significant smooth- 
ing used to generate the temperature map. The overall structure is 
strikingly similar to the predictions of the simple model considered 
in the previous section. It is therefore reasonable to assume that 
this cool gas is the low entropy gas transported from the core of the 
cloud to the interface where it expanded and cooled to the observed 
low temperature. 



4 DISCUSSION AND CONCLUSION 

Cold fronts were also identified in the full (cold dark m atter 
+ fluid d ynamics) simulations of iNagai & Kravtsov (2003) and 
iBialek et'al. 1 2002,) . Compared to these simulations our model is 
intentionally kept simple, using rather idealised initial conditions. 
This simplicity, however, allows us to see the evolution of the cloud 
characteristics in a physically more transparent and controlled sit- 
uation. The results sh own above are in broad a gree ment with the 
sugges tions made by iNagai & Kravtsov! i2003h and iBialek et alJ 
i2002h . although the detailed description of the cloud evolution is 
different. As such, our simulations are not meant to replace full 
blown 3D cluster merger simulations, rather, we use them as a 
tool to point an important mechanism at work inside ram pressure 
stripped sub-clusters, which, while present in larger scale simula- 
tions, had not been identified in previous publications on the sub- 
ject. 

We show that slow gas motion inside a gaseous subcluster 
moving through a uniform medium transports gas from the central 
parts of the cloud towards the contact discontinuity, which sepa- 
rates the subcluster gas from the external gas. In a typical cluster 
or group of galaxies, the entropy of the gas in the central region 




Figure 2. Projected temperature distribution (in keV). The contours are the 
isophotes of the X-ray surface brightness in the 0.5-3 keV energy band. 
The outermost contour approximately traces the discontinuity, discussed in 
IVikhlinin et alJ The extension of the cool region outside the dis- 

continuity (outermost contour) is due to the significant smoothing of the 
temperature map. The image is 26 by 26 arcminutes. 



is significantly lower than the entropy of the outer layers. The gas 
transported from the core region will undergo adiabatic expansion 
and cool to temperatures less than its initial temperature. This cool 
gas may form a cool and bright tip at the leading edge of a mov- 
ing cloud. Such a situation can be characteristic of a subcluster or a 
group that still retains gas in its gravitational potential, but spends 
sufficient time moving through the cluster gas so that motions in- 
side the cloud are fully developed. 

Assuming that the original subcluster had a metallicity gradi- 
ent with the most metal rich gas near the center, the transport of the 
core gas towards the interface should also enhance the abundance 
contrast across the discontinuity. 

As mentioned before, the simulations presented in f|2|ideal- 
lized for the purpose of isolating the process we were interested 
in studying in this paper. Many more levels of complexity can be 
added to tailor the simulations to specific cluster merger scenarios. 
For example, it is easily possible to simulate such a merger by run- 
ning two beta model atmospheres into each other. In this way, one 
can capture the additional effects of varying ram pressure due to 
the stratified cluster atmosphere and of tidal effects. 

To verify that the effects discussed above act in such a more 
general scenario, we ran several simulations of two beta model 
atmospheres colliding, with one cluster being significantly larger, 
more massive, and hotter than the other. Indeed, the same mecha- 
nism discussed above operates in such cases as well. In addition, 
more complex dynamical effect occur in the stripped wake, es- 
pecially immediately after the passage of the smaller cluster core 
through the bigger one. This interaction also induces some in- 
ternal waves which produce interesting observable features inside 
both atmospheres (comparable to the effects of sloshing mentioned 
above). Additionally, adiabatic expansion of the material as it ex- 
its the larger cluster after passage leads to a significant cooling of 
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the material, enhancing the effect described above even more: The 
lowest entropy and thus coolest material lies at the front of the sub- 
cluster. 

Such simulations introduce a number of additional free pa- 
rameters, and it is beyond the scope of this paper to map out all 
of this parameter space, in particular since even such a simulation 
will be artificial in that it neglects the evolution of the underlying 
dark matter, truncation effects of the cluster atmospheres, and non- 
sphericity of the gas clouds. We note, furthermore, that reproducing 
a specific observational scenario in detail, such as the example of 
Abell 3667, is a non-trivial exercise in fine tuning. For example, 
the estimated low Mach num ber of around M ~ 1 seen in this 
object jVikhUnin et alj|2001ij) is difficult to reproduce by simple 
free-fall simulations of two clusters, which ultimately lead to Mach 
numbers of around M 3 near the impact of the two cores, corre- 
sponding to the typical depth of the potential well of a King model 
I Kina 1962). The frequency of incidence of such low Mach number 
merger shocks is a question that can only be answered by detailed 
studies in the context of structure formation simulations. 

Thus, many other effects will play a role in determining the 
observational appearance of subclusters in ram pressure dominated 
environments. For example, a subcluster which has already passed 
through the dense core of the main cluster will have been stripped 
of most of its outer gas. Upon leaving the core on the other side of 
the main cluster, the subcluster will propagate into a lower pressure 
and lower density environment, thus expanding adiabatically and 
cooling to form a cold cloud, as observed in many cluster environ- 
ments. Mergers may also induce sloshing of low entropy gas in the 
potential wells of the cluster cores and cause the displaceme nt of 
the cool gas relative to the equilib rium position I Markevitch et alj 
I2OOI '. 2002:'C hurazov et all2003l) . Nonetheless, the processes de- 
scribed in this paper will operate under most of these conditions 
as well, leading to the transport of cold material to the front of the 
cloud on time scales comparable to the lifetime of the subcluster. 
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tion and NASA contract NAS8-39073. 
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